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Abstract 

We present the effective lagrangian for low energy and momentum spin waves in 
canted phases at next to leading order in the derivative expansion. The symmetry 
breaking pattern SU{2) ^ 1 of the internal spin group and that of the crystallographic 
space group imply that there is one ferromagnetic and one antiferromagnetic spin wave. 
The interaction of the spin waves with the charge carriers is also discussed for canted, 
ferromagnetic and antiferromagnetic phases. All this together allows us to write the 
doping dependence of the dispersion relation parameters for doped manganites. We 
point out that the spin waves posses distinctive characteristics which may allow us 
to experimentally differentiate canted phases from phase separation regions in doped 
manganites. 
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1 Introduction 



Canted phases are magnetically ordered states with non-collinear magnetisations. These 
configurations appear in quantum Hall double-layer systems and in the conducting regime 
of double exchange models 0, |, || , where the local magnetisations arrange in two sublattices 
with magnetizations pointing to different (but not opposite) directions. The double exchange 
models are believed to provide a good description of doped manganites |^ , which are receiv- 
ing quite a lot of attention lately p, 0, ||, |1^, [n|, |12|, [1^. Doped manganites present a 
non-trivial interplay between their magnetic and conducting properties |15|, |l^, [l^, |18|, |19 



which leads to a rich phase diagram. The transitions between the different phases in terms 
of the doping have been extensively studied . 



The most studied transition is that from an antiferromagnetic insulating phase, at zero 
doping, to a ferromagnetic conducting phase as the doping grows It is not clear yet 

if the region for intermediate values of doping corresponds to a canted phase or to a phase 
separation region | 



23| , |2^ , |25| . It is our claim that the study of the spin waves in such 



materials may shed light to this question. Since the spin waves are low energy excitations in 
a magnetically ordered material, they are sensible to the main features of the phase diagram. 
The spin waves have indeed been studied recently in connection with these materials p6|, 

The low energy and momentum dynamics of the spin waves is so much constrained by 
the symmetries of the system that a model independent description is possible in terms of a 
few unknown parameters. For canted phases, the spontaneous symmetry breaking pattern 
of the (internal) spin symmetry is SU{2) 1, instead of SU{2) U{1) like in ferromagnets 
or antiferromagnets. (Strictly speaking the symmetry breaking pattern is SU{2) —>■ Z2, the 
center of the group. However, since we will not be concerned with global properties, neither 
of the group nor of the coset manifolds, using 1 instead of Z2 does not modify our discussion 
at all). As a consequence of the Goldstone's theorem there will be gapless excitations in 
the spectrum (Goldstone modes) [0, the so called spin waves. A very efficient way to 
encode the spin waves dynamics is by using effective lagrangians. 

Effective lagrangians for Goldstone modes are known since the late 1960s and they 
have been extensively used in pion physics during the last decade It was suggested in 
ref. ||3^ that they may also be useful in Condensed Matter systems. A detailed construction 
of the effective lagrangians for ferromagnetic and antiferromagnetic spin waves has already 

40| for a recent review and |^T], ^ for non-trivial applications). 



been presented in (see 

A general formalism for the spin waves in canted phases is presented in Section 2, where 
we construct an effective lagrangian at next to leading order. An intuitive separation of the 
spin wave field in one ferromagnetic and one antiferromagnetic component is also presented. 
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Since canted phases appear in the conducting regime of doped manganites, the couphng of 
spin waves to charge carriers is relevant. This is discussed in Section 3, where we obtain 
an effective lagrangian for this couphng in the three different phases: canted, ferromagnetic 
and antiferromagnetic. In Section 4 we use the previous resuhs for the different phases 
of doped manganites in order to obtain the exphcit dependence of the dispersion relation 
parameters on the doping, which is given in formulas ( ^.191 ), ( [4 .211 ) and ( ^.221 ) for the canted, 
ferromagnetic and antiferromagnetic phases respectively. In Section 5 we present a plot of 
the doping dependence of the velocity and the mass of the spin waves for the different phases. 
We also explain in Section 5 how our results on spin waves can be used to experimentally 
disentangle canted phases from phase separation regions. We summarize our conclusions in 
the last section. Some properties and calculations, related to a loop integral, are relegated 
to the Appendix in order to keep our arguments clear. 

In order to simplify the notation we take h = c = 1, which leads to a relativistic notation. 
Hence, we use x = (t, x), g = {u, k) and subindices = 0, 1, 2, 3, where the zero stands for 
the time component. Space indices are denoted by i = 1, 2, 3. 



2 Effective Lagrangian for Canted Phases 

In a previous paper we obtained the phase diagram for doped manganites, where a 
rich set of magnetically ordered phases appeared. The magnetically ordered configurations 
break spontaneously the SU{2) symmetry of the theory (the continuum double exchange 
model) down to the ground state symmetry, U{1) for the ferromagnetic and antiferromagnetic 
configuration, and 1 for the canted configuration, because of the non-collinear character of 
the latter. In this situation the lower excitations of the system are the spin waves, which 
turn out to be the Goldstone's modes associated to the spontaneous symmetry breaking in 
magnetic systems. 

We have already carried out an extensive study for the ferromagnetic and antiferromag- 



netic spin waves in crystalline solids in |39]. This formalism assumes the existence of a gap 
in the excitation spectrum, which permits the construction of an effective lagrangian for the 
spin waves as an expansion of local terms suppressed by the gap. In the case that there are 
additional degrees of freedom with energy smaller than the gap, they should also be included 
in the effective lagrangian. This is the case of charged carriers in doped manganites, which 
we will discuss in section 3. In this section, we restrict ourselves to the generalization of the 
formalism presented in [R^ to the case of spin waves in canted phases. 
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2.1 Effective fields and symmetries 

The (internal) spin symmetry breaking pattern, SU{2) —>■ 1 for the non-colhnear canted 
configurations, determines that the basic field which represents the Goldstone modes (spin 
waves) may be chosen as a matrix V{x) G SU{2)/1 = SU{2) [^. After determining the 
transformations of this field under the symmetry group of the system we can build an effective 
lagrangian from which the spin wave dynamics can be derived. The transformations under 
the SU (2) spin symmetry read 

Vix) ^ gV{x) , geSUi2). (2.1) 

The transformations under the crystallographic space group reduce in the continuum to 
the primitive translations and the point group. Since the local magnetizations in the canted 
phase point to two different directions depending on the site, both of these symmetries are 
broken by the ground state. This must be refiected in the transformation properties of V{x). 
For definiteness, we shall take the local magnetizations in the 1 — 3 plane in the spin space, 
in such a way that the even and odd lattice magnetizations form an angle of 6/2 and —6/2 
with the 3-axis respectively, and can be mapped into each other by a rotation of vr around 
the 3-axis. 

Ml = M(sin(0/2),O,cos(e/2)) 

M2 = M(-sin(^/2),O,cos(0/2)) (2.2) 

For simplicity, we shall also assume a primitive cubic lattice {PmSm) although the ana- 
lysis can be carried out in a similar way for any crystallographic space group. The point 
group mSm is generated by the transformations C22, C2y, C2a (a 2- fold axis in the direction 
(1,1,0)), C^i (a 3-fold axis in the direction (1,1,1)) and the spatial inversion /. These 
transformations can be separated in two groups, on the one hand {C2z,C2y,C^i}, which 
transform points inside each sublattice, and on the other {C2a,I}, which as the primitive 
translations, r, transform points from the even sublattice to the odd one and vice-versa. 
Thus the transformation of the spin wave field V{x) under this group is given by: 

^■.{C2.,C2y,C^i}: Vix) ^ g^Vix) 

^■.{C2a,I}: Vix) ^ g^Vix)R , R = e-'-'\ (2.3) 
r: Vix) Vix)R 

where g^ is the SU (2) transformation associated to the point group transformation and R is 
a matrix which interchanges the magnetisation between sublattices. 

Notice that by combining the transformation of the field Vix) under {C2a, 1} with the 
translations in {rC2a, tI} we can eliminate the additional factor R in those point group 
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transformations. Since, in addition, the factor can be re-absorbed by a SU{2) transfor- 
mation we only have to care about the transformations of the derivatives as far as the point 
group is concerned. 

Finally, under time reversal V{x) transforms as 

T: V{x) V{x)C , C = e-'^'^\ (2.4) 

We are now in a position to construct the effective lagrangian order by order in derivatives. 
In order to do that we consider the following element of the Lie algebra of SU{2) |3^ : 

V\x)td^Vix) = b-^{x)S+ + b+{x)S-+bl{x)S\ (2.5) 

This term, and consequently the coefficients &Ji(x), are invariant under the SU{2) transfor- 
mations ( p.l|) . Under the point group we only need to consider the transformation of the 
derivatives in 6^ (a;), which correspond to the space-time indices fi. The transformations 
under primitive translations are given by 

, Z 7 <^-^> 

and under time reversal 

( S " "5 (2-7) 
where tfi stands for the transformation of the index n under time reversal T. 

2.2 Relation with ferro and antiferromagnetic spin waves 

Before writing down the effective lagrangian, let us discuss a suitable decomposition of 
V{x) which illuminates the relation between canted spin waves and the usual ferromagnetic 
and antiferromagnetic ones. In the way we have chosen the direction of the magnetizations 
in each sublattice ( |2.2|) , it is clear that the projection on the third direction is ferromagnetic, 
whereas the projection on the 1 — 2 plane is antiferromagnetic. This suggests that we may 
separate the spin wave field into components perpendicular to the third axis and to the plane 
1 — 2 respectively. Group theory allows us to implement this easily. Indeed, an element of 
the group, V{x) G SU{2), admits a unique decomposition in terms of an element of a coset, 
U{x) G SU{2)/U{1), and an element of the corresponding subgroup, H{x) G U{1), such that 
V{x) = U{x)H{x), with 

U{x) = exp\^\7T-{x)S+ + TT+{x)sj\ G SU{2)/U{1) 



H{x) = expl'-^7c\x)S^\ G f/(l). 



(2i 
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where S± = ± iS^ and are the SU{2) generators, 7r^(x) = (tt ^{x)±m^{x)) 1^/2 and 
7r^(x) are the spin waves fields, and and /s are dimensionful parameters representing the 



spin stiffness. This imphes that the element of the Lie algebra in ( |2.5| ) can be written as 

V\x)id^V{x) = H\x){U\x)id^U{x))H{x) + H\x)idf,H{x). (2.9) 
Upon using for U'^ {x)id ^JJ {x) a similar expression to that in ( ^.5] ) we have, 

U\x)id^U{x) = a^{x)S+ + a+{x)S^ + al{x)S\ (2.10) 
This decomposition translates to the coefficients &J^(a;) in (|2.5|) as follows: 

b-^{x) = e-^^^-'^^^lf^a-{x) 

h+{x) = e'^"'(")/^3^+(x) (2.11) 
hl{x) = alix)~V2dyix)/fs. 



Recall finally that the expansion of UHdnll in spin wave fields reads 



J n 



(2.12) 



2.3 Effective lagrangian 



In order to construct the effective lagrangian let us begin by considering terms with time 
derivatives. It is then clear that we can build a term with a single time derivative, 

bl = al-^do7r'' r^al (2.13) 

which contributes to the dynamics of ir^^x). Nevertheless since this term only contains a 
total derivative on tt^{x) the first contribution to the dynamics of this field comes from 



~ -j^do-K^doTT^ (2.14) 

/3 



where we have made explicit the quadratic term in tc^Ix). 

Regarding the spatial derivatives there are no invariant terms with a single spatial deriva- 
tive. Then the first invariant terms have two space derivatives, and they read 



KK + KK (2.15) 
2 

p2" 
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where again we have made exphcit the quadratic dependence on 7t^{x) in the last term. 

Unhke the terms with time derivatives, the terms with spatial derivatives produce a 
leading order contributions for 7c^{x) and n^lx) at the same order. Let us call it 0{p'^). 
Equations (|2.13|) , ( p. 14 ) and ( 2.15 ) provide the dispersion relations for the spin waves, which 
indicate how time derivatives must be counted with respect to space derivatives. Namely, 
a time derivative on tt^{x) must be counted as O(p^), whereas a time derivative on 7r^(a;) 



must be counted as 0{p). This implies that the term 6q 



v^9o7rV/3~0(p2) + 0(p), 



i.e., it contains terms of first and second order, which must be taken into account in the 
construction of the effective lagrangian. This is in fact a remarkable difference with respect 
to the ferromagnetic and antiferromagnetic case, where each invariant term has a unique 
size. 

Then, putting together all the terms above, the most general effective lagrangian at order 
O(p^) we can construct, with the standard normalizations, reads 



C(x) 



f 

J 77 



2 ° 



Bhihl 



KK 
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(2.16) 



If we expand it up to three fields, it reads 



Cix) 



C 

— {di7r+di7r+ + diTT~diTT' 



1 
2 

iV2C 



V 

~2 



(2.17) 



iV2v'f3 



3 ^/tt 

The first two lines correspond to quadratic terms in the fields, which yield the free 
propagation of the spin waves and give the dispersion relation for each of them. Whereas 
these terms lead directly to a wave equation (Klein-Gordon type) for tt^^x), as expected for 
an antiferromagnetic spin wave, the equation for 7t^{x) turns out to be non-diagonal. The off- 
diagonal terms are due to the existence in the effective lagrangian of the term {bfbf + h^h^). 
In the ferromagnetic and antiferromagnetic cases this term does not appear because the 
unbroken U{1) subgroup prevents it. In order to diagonalize the quadratic 7r^(x) terms, we 
perform the following Bogolyubov transformation: 



TT X 



Im'B 



7r+(x) 



'm'B 



1 



TT (x) 



1 



2m' 



- C^. (2.18) 



In terms of the new variables the lagrangian ( |2.17|) reads 

1 _ 1 v'^ . 

V J u * * 2 2 



iV2C 



dill djir )7r^ 



7^+9i7^-)9i7^^(2.19) 
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which yields a Schrodinger equation with a mass m! for the new field 7r+(x). Therefore, 
as it was expected from the decomposition made in (|2.8| )- (|2.11| ) the field 7r+(x) describes 
one ferromagnetic spin wave, with a quadratic dispersion relation, and 7r'^(x) describes one 
antiferromagnetic spin wave, with a linear dispersion relation. 

This result is in agreement with previous theoretical [IH and recent experimental |^| 
works, and in line with the general counting of Goldstone modes in non-relativistic systems 
stated in |^ (see also [^). The general statement is that there exist as many real fields 
representing the Goldstone modes as broken directions in the symmetry group (three in 
our case, because of the non-collinear nature of the canted configuration). The space-time 
transformations for these fields determine if they verify a wave (Klein-Gordon, leading to a 
linear dispersion relation) or a Schodinger (quadratic) equation of motion, with the constraint 
that in the case of a Schrodinger equation a complex field, and therefore two real ones, is 
necessary to represent a single physical mode (the two real fields behave like canonical 
conjugate degrees of freedom). With this argument in mind for the canted spin waves we 
can only get either three linear branches or one linear and one quadratic branches, which 
turns out to be the correct answer in our case. 

At next to leading order, O(p^), besides those terms coming from 6q6q already considered 
in ( |2.16| ), we find the following terms: 

blibtbt + brbr) 

blb-bf (2.20) 

%"i "i ■ 

The term i{b^dobf — bfdob^) ~ blb^bf at this order, since at O(p^) we only have to consider 
time derivative acting on 7r'^(x). However at higher orders the terms obtained from those 
invariants are different from each other. 



2.3.1 Coupling to a magnetic field 

The most important source of magnetic coupling in a spin system is the Pauli term, the 
introduction of which in the effective theory was extensibly discussed in |^ . The outcome 
is that the Pauli term can be introduced by just replacing the time derivative by a covariant 
derivative in the following way: 

do ^ Do = do- ifimSB. (2.21) 
After introducing the covariant derivative the equation (|2.9| ), for time derivatives, reads 
V^{x)iDoV{x) = H\x){U\x)iDoU{x))H{x) + H\x)idoH{x). (2.22) 
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Thus, after introducing the magnetic field the effective lagrangian is constructed with 
the expressions such that the magnetic field only modifies a^{x) and clq^x), given by 



UUDoU = - — 

J TT 



+ 

+ 



UdoTT^ - fir, 



1 



-TT 



S. 



^4 



(2.23) 



X] 



The time derivative on vr^(x), as well as the terms with spatial derivatives hf{x) and 6f( 
remain unchanged by the presence of the magnetic field. 

It is very easy to see that, at the lowest order, the dispersion relation of the antiferro- 
magnetic branch, given by vr^(x), is not modified by the introduction of a small magnetic 
field in any direction, in particular in the direction of the staggered magnetization. 



2.3.2 Ferromagnetic limit 



In the ferromagnetic limit the local magnetizations are pointing in the third direction all 
over the crystal. Hence, an unbroken f/(l) symmetry remains, and the spin waves are repre- 
sented by a field belonging to the coset SU{2)/U{1). This field can be easily obtained from 
the decomposition (|2.9| ) of the canted case by taking H{x) = 1, or, equivalently, tt^{x) = 
in Hence, V{x) simply reduces to U{x). Furthermore, because of the remaining 

U{1) symmetry terms like {afaf + a^a^) are forbidden, and hence the quadratic part of the 
lagrangian does not contain off-diagonal terms. Therefore the Bogolyubov transformation 
( p.l8| ) is not necessary anymore. 

In terms of a^(a;) the effective lagrangian for the ferromagnetic spin waves reads 



C(x) 



J TT 



1 



1 3 

-an 

2 ° 2m' 



a, a] 



And after expanding it in terms of spin wave fields in (|2.12| ), 

1 



Cix) 



2m 



-diir di 



TC 



(2.24) 



(2.25) 



which corresponds to one spin wave with a quadratic dispersion relation. 



2.3.3 Antiferromagnetic limit 

In the antiferromagnetic limit the local magnetizations are pointing at opposite directions 
in each sublattice along the first axis (5*^). As for the ferromagnetic case an unbroken 
U{1) symmetry remains, and the spin waves are represented by an element of the coset 
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SU{2)/U{1). In order to simplify the computation we will rotate the internal space reference 
frame in such a way that the third direction, instead of the first, lies along the staggered 
magnetization direction (we perform the rotation 1 ^ 3 — > 2 in all the indices). With this 
choice the spin wave field is determined from ( p.9|) by setting H{x) = 1, or, equivalently, 
7r^(x) = in (|2.11|) . V{x) reduces to U{x), and the remaining U{1) sjTumetry prevents 
the non-diagonal terms, {afaf + a^a~), from appearing like in the ferromagnetic case. In 
addition to that, now C acts as the matrix which interchanges the magnetizations between 
sublattices, which forbids the term with a single time derivative, al{x), to appear in the 
effective lagrangian. 

The effective lagrangian for the antiferromagnetic spin waves is given by 



And after expanding it in terms of spin wave fields in ( |2.12| ), 



(2.26) 



(2.27) 



which describes two spin waves with a linear dispersion relation. These two branches are 
splitted by the introduction of a small magnetic field in the third direction, the direction of 



the staggered magnetization [R^ 



3 Interaction with charge carriers 

Canted phases are known to support conductivity. Then it is important to elucidate 
which kind of interaction mediates between the spin waves and the charge carriers. In order 
to address this question in a model independent way, we would need an effective field theory 
description of the latter. However, to our knowledge, there are no general rules on how to 
build such an effective theory, which may depend strongly on the particular material we wish 
to study. We shall then restrict ourselves to present an effective theory based on a particular 
model which successfully describes canted phases and has applications to doped manganites, 
the continuum double exchange model p2| . 

At first sight one may think of describing the charge carriers by an effective fermion field 
which varies slowly through the material and couples to the local magnetization. However, in 
a canted phase the local magnetization changes abruptly from the even to the odd sublattice, 
which means that we shall need two magnetization fields Mi(x) and M2(x) in the even and 
odd sublattices, and hence a single slowly varying fermion field is not enough to have a 
consistent description. We need at least two slowly varying fermion fields ^Jii^x) and ip2{.x), 
coupled to the magnetization in the even and odd sublattices respectively. 
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The interaction lagrangian of the model reads 



q2 ^ 

(1 + it)idQ + ^ + + J//-Mi(x) 



+ t (■?/'i(x)-?/'2(x) + iIJ2{x)iIJ2{x 



a- 

1 + ie)idQ + ^ + + JH^yi-2{x 



%l)2{x) 



(3.i: 



where t corresponds to the amphtude of probabihty that the fermion changes the sublattice 
and Jh is the Hund couphng between the fermion fields '4)i{x) and ip2{x) and the magnetic 
moment in each sublattice Mi(x) and M2(a:;) respectively. An estimation of our parameters 
is given by t ~ zt^, Jh ~ Jh 2m ~ where a is the lattice spacing, z = 6 is the 

coordination number and the superscript / means the analogous lattice quantity. In order 
to have conduction when t ^ only, the chemical potential fi is required to lie below the 
lowest energy of the band for t = 0. 

The spin waves are fluctuations of the magnetically ordered ground state, and they are 
included in the previous fields. We can separate the contribution of the spin waves from 



that of the ground state Mi and M2 by writing Mf (x) = R^{x)M^ 



1,2), such that 



the matrix R^{x) corresponds to the spin wave fluctuation in the adjoint representation of 
SU{2). Using the scalar product properties the interaction term can be written 



-M,(x) = yi?^(a;)M^ = V{x)^V\x)Mi 



1,2), 



(3.2) 



where the matrix V{x) represents the spin waves in the fundamental representation of SU{2). 
This expression suggests the following change of variables for the fermionic fields 



^lj,{x) V{x)iJi{x) (z = l,2). 

In terms of the new fermionic fields the lagrangian (|3.1| ) reads 

^l(x) iplix) 



Cix) 



O^' + 




(3.3) 



(3.4) 



where 
09' -- 



;i + ie)idQ + d^/2m + /i + ^o-Mi 



t 



[I + ie)ido + 52/2m + /i + ^o-M2 



(3.5) 



is the contribution of the ground state, and 








= (vUdoV) - — [{td,, (VhdiV)} + iV\d,V){V\diV) 



(3.6) 
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contains the interaction with the spin waves. The curly brackets { , } stand for the anti- 
commutator. Taking into account the decomposition (p.5|) the operator O^^ can be expressed 



in terms of the fields as follows: 



2m 



2m 



{id,M}]S- 



(3.7) 



A compelling expression for the coupling of the spin waves can be written by noting that 
the expressions ( |3.4|) -( pr7D are equivalent to the introduction of a covariant derivative, 



id,, 



iD^ = id^ + {V^id^V) = id^ + h-s+ + h+s_+hls^ 



{3.t 



in ( |3.5|) and dropping O^'^. 

Since the spin waves are fluctuations of long wavelength and the interaction with the 
fermionic fields contains derivatives, this interaction will be small. In this situation the 
problem is reduced to calculate the interaction of the spin waves with the eigenstates of 
( p.5|) perturbatively. The four eigenstates can be obtained by considering the following 
change of variables: 



X2{x) 



pt 



J_ q + Q^ + Q q-Q-Q 
V2 \ q-Q^ + Q -q-Q + Q 



(3.9) 



where q is an scalar parameter and and Q = are matrices in the Lie algebra of SU (2) 
given by 



\ 



e+ + 7 + cos 2 
2^1 



\ 



e_ + 7 - cos 2 1 7<i 1 



2e_ 



Q 

with 



sm 



2e_ (^e_ + 7 - 



sm ■ 



cos ; 



2e+ (e+ + 7 + cos ■ 



e 

^ , cos - + sm - 
2 \ 4 4 



D_ — > cos - D+ + sm - D_ 



\ 



e+ + 7 + cos ■ 
2^1 



2e_ 





( ^ 




cos- 











e± = W 1 + 7^ ± 27 cos ■ 



e 



2t 



1 



(3.10) 



(3.11) 



2 ' ' \Jh\M' 
M = I Ml I = IM2I =3/2, and 6 the angle formed by the ground state magnetizations Mi 
and M2. 
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After the change of variables the lagrangian, written in terms of the new fields Xi{^) 
X2{x), reads 



C(x) 



Xi{x) xlix) 



Li 



+ 



W21 W22 




L2 

where the interaction with the ground state is diagonal, and it is given by 

\Jh\m I ~ e 

1 / 1 + 7"^ ± 27 cos - 



(3.12) 



L2 = (l + ^e).5o + ^ + /.-^ 





1 + 7^ =F 27 cos -. 



(3.13a) 
(3.13b) 



The interaction with the spin waves, given hj W = PO^^P\ reads 



W^2 

W22 



Ml = -g[o^"',g]-(go^"'g + Qo^"Q) 



(3.14a) 
(3.14b) 
(3.14c) 



where the square brackets [ , ] stand for the commutator and O^'^ is given in ( p.7| ). 

In the relevant materials that we have in mind, the hopping amplitude, t, is much smaller 
than the Hund coupling, Ju-, i-e. 7 -C 1. In this case, the two higher states of (|3.5| ), denoted 
by X2(a;), lie far away from the two lower ones. In fact, the ratio of energies is of order 7. If 
we are only interested in transition energies ~ t we can safely integrate out the states X2(a;), 
obtaining the following lagrangian for the two lowest states, X\ip^\- 



^eff = x\{x) {Li + Wu) Xi{x) - x\{x)Wu 



L2 + W, 



-W2iXi{x). 



(3.15) 



22 



The second term is of order 7 with respect to the first one. Indeed, we consider low 
incoming energy and momentum with respect to the two lowest states, namely, 
Li ~ (1 + ie)ido + df /2m + /i + \ Jh\M/2 ~ t, and the spin wave interaction Wij ~ t. Thus 
1/1/2 ~ l/t, such that the second term in ( p.l5|) is order •yt, which means that the field X2{x) 
decouples, and the effective lagrangian reduces to 



x{{x) (Li + Wn + 0{jt))xi{x). 



(3.16) 



In order to complete the effective lagrangian we must consider the leading order in 7 for Wu 
in ( |3.14| a), which corresponds to take the right limit in (|3.10| ) for the parameter q and the 
matrices and Q. 
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Therefore, the effective lagrangian is given by (|3.16|) , where the interaction terms, which 
come from (|3.14| a), (p.7|) and ( p.lO|) , are given by 
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1 ^ /. 1 



hthT + ) + ^ cos ^ f - —{id,, hi} 



1 
2m 



+ - sin - [h-, - —{id., h-} + 6+ - 6+} ) (5+ + 5_) . (3.17) 



Upon expanding it up to two fields, using ( p.l2| ) and (|2.11|) , we finally obtain 
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2m/2 



J TT 



1 e 



2/1 



+ 



fs V 2m 



1 e 

H sm - 

2/./3 2 



(3.18) 



2m 



+ ^9i7r^(9i7r - diTT^ 
2m 



3.1 Coupling to ferromagnetic spin waves 

The interaction of charge carriers with ferromagnetic spin waves can be considered as 
a limit of the canted configuration. In order to do that we must take the limit ^ in 
( p.lOj ), which yields a very simple expression, independent of 7, for the parameter q and the 
matrices and Q which determine the change of variables ( |3.9| ), namely, 

q = I = Q = S'. (3.19) 

Following the considerations in subsection 2.3.1, we also must take H{x) = 1, or equivalently 
n^{x) = 0. 

Because of the remaining unbroken symmetry the SU{2) transformations on the spin 
wave fields are realized by a non-linear U{l)iocai gauge group, which allows us to write the 
lagrangian in a manifestly gauge invariant way 

\Jh\M 

I UiU^ — a - a; \ + 11, + ■ 

2m 



Li + Wii = iDo + ^ (DiDi - a^af) + /i + + 2tS^ 



(3.20) 



where iD^^ = id^ + a?^/2. Notice that this implies that some of the couplings are fixed by 
the symmetry. In fact only the coupling a^af is model dependent. This is analogous to 
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what happens in the pion-nucleon lagrangian where one of the couphngs is fixed by chiral 
symmetry Recall that the transformation properties under U{l)iocai are the following: 



^'-^ e^*^(")aj(x) (3.21) 



Finally, in terms of the spin wave fields and up to two fields W\\ reads 
2mf4 



1 

+ 



(3.22) 



which corresponds to the limit ^ of the canted expression (| 3T 

3.2 Coupling to antiferromagnetic spin waves 

In the antiferromagnetic case we have to consider two situations. The first one corre- 
sponds to the insulating phase, where there are no charge carriers to couple with. The 
second situation corresponds to the antiferromagnetic conducting phase. We shall describe 
this second situation below. 

The interaction of the charge carriers with antiferromagnetic spin waves is given by the 
canted case in the limit 6 —>■ n. According to the discussion in subsection 2.3.2 this limit 
is a little bit more involved than for the ferromagnetic case, since we also must rotate the 
reference system (1 — 3 — 2). Then the expressions for the parameter q and the matrices 

and Q, which determine the change of variables ( |3.9D , read 



gt = ' ^ 2^ V25^ (3.23) 

V 1 + 7^ + 7Vl + Y 

Q = 0, 

where the right limit gives the leading dependence on 7 which we shall use to calculate the 
interaction. 

Similarly to the ferromagnetic case we must take H{x) = 1, or 7[^{x) = 0, and therefore 
the remaining unbroken symmetry determines the gauge invariance structure for the effective 
lagrangian, given by 



Li + Wn =tDo + — [D,D, - a^a+j + /i + (3.24) 
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where now iD^ = id^ + a^S^. Notice again that this imphes that some of the couphngs are 
fixed by the symmetry. In fact only the couphng a^af is model dependent. This is analogous 
to what happens in the ferromagnetic case discussed before and, hence, also analogous to 



the case of the pion-nucleon lagrangian ^J] . Recall that the transformation properties under 
U{1) local are now the following: 



aj(x) 



(3.25) 



Finally, in terms of the spin wave fields and up to two fields Wu reads 
1 



2m/2 



84 n (9, 77^ 



1 

J TT 



TT 



5^(3.26) 



which corresponds to the limit ^ tt of the canted expression ( |3.18| ). However because of 
the rotation of the reference system it is easier to obtain ( p.26| ) as a limit of ( |3.17|) rather 
than (|3l8| ), because the direction of the symmetry breaking has not been taken explicitly 
yet. 



4 Spin Waves in Doped Manganites 

In the section 2 we developed a general formalism which provides the effective lagrangian 
for the spin waves generated by any model in a canted, ferromagnetic and antiferromagnetic 
ground state as long as the model is invariant under SU (2) transformations. In this section we 
are going first to particularize this effective lagrangian to the case of the spin waves in doped 
manganites. Next, we will include the interaction of spin waves to charge carriers worked 
out in section 3 and calculate the doping dependence of dispersion relation parameters. 



4.1 Spin waves from the Heisenberg hamiltonian 

In the double exchange models the interaction between the core spins in the t2g-bands of 
the manganese atom is described by an antiferromagnetic hamiltonian. Since the value of 
the core spins is 3/2 their motion is slow and can be approximated by classical magnetization 
fields on the lattice. Furthermore, for the low energy and momentum region the lattice fields 



can be further approximated by continuum fields. In p2| we considered a static Heisenberg- 
like interaction, which only provides the relevant contribution to the ground state energy. 
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Here we shall introduce a derivative expansion of the Heisenberg hamiltonian, which also 
takes care of the low energy and momentum excitations. These derivative terms in the 
second quantization language read 



H 



d^^^^dMi{x)dM2{x), 



(4.1) 



where Jaf ~ ^Jaf/^^ the superscript / represents the lattice Heisenberg coupling. The 
local magnetizations for each sublattice, Mi(x) and M2(x), are given as fluctuations of the 
ground state configuration, as mentioned in section 2, that we can write ipi{x) = V{x)ipf'\ 



M,{x) = (^f^Vt(x)S\/(x)(^f ^ = tr (v\x)SV{x)Pi) = Rl{x)M^ 



(4.2) 



where Pi is a projector in the direction of the ground state magnetization in each sublattice 
z = 1,2. 

This hamiltonian only generates terms with spatial derivatives in the spin wave's effective 
lagrangian. In order to introduce the temporal term let us consider it written in terms of 
the total, S(x), and staggered, r2(x), magnetizations 



H = - 



'^^^''^■dii:{x)di'E{x) - :l^f^din{x)din{x) 



2z 



8z 



(4.3) 



where 



S(x) 



[Mi(x) +M2(x)) 



n{x) = Mi(x) - M2(a;) 



Mcos' 



Q = 2Msin-. 

2 



(4.4) 



This corresponds to the spatial derivatives terms in the effective lagrangian for ferromag- 
netic spin waves in terms of the total magnetization, S(x), and for antiferromagnetic spin 
waves in terms of the staggered magnetization, fi{x), |^5|. Following this identification we 
shall choose as temporal terms those which complete these hamiltonians. Then the effective 
lagrangian from the Heisenberg contribution reads 



X 



a 
+ 



i- rfAS(a;, A) ido^{x, A) x dx^ix, A)) + ^^^di^{x)d,^ix) 

1 o 7-2 802 

-9on(x)9of^(x) - -^l-^^,nix)^^nix) 



UJafo^^'^ 



(4.5) 



where S(x, A) is an extension of the total magnetization field which verifies S(x, 0) = S 
and S(a;, 1) = S(a;). This is equivalent to introduce an extension 7r°(a;, A) for the spin 
wave fields such that 7r'*(a;, 0) = and tt°'{x, 1) = 7t°'{x). A simple extension valid for our 
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purposes is 7r"(x, A) = Xtt°'{x), which allows us to write the effective lagrangian for the canted 
configuration in terms of the spin wave representation used in section 2, 



x) 



2E 



+ 



A 



3J^jra^r2^ 3 3 



■ ibtbt + brb;) 



(4.6) 



We have dropped terms with two time derivatives acting on vr (x) since they are sub- 
leading in the canted and ferromagnetic phases. However, they are not so in the antiferro- 
magnetic phase and will have to be restored in order to take the antiferromagnetic limit. 
Notice that ( [4.6|) provides particular values for the constants B, C, /| and v"^ in the 
general formula ( p. 16 ). 



4.1.1 Ferromagnetic configuration 

In the case we have the ferromagnetic configuration the limit is taken very easily, since 
in this case E ^ M and ^ 0. Since the time evolution is already described by the term 
with a single time derivative we can drop the two time derivatives term in ([4.6|) , which yields 
the effective lagrangian 



£«( 



2M 



x) 



1 3 

2^0 + 



JApa^M 



-a, a 



(4.7) 



Notice that the mass term in ( [4.7|) has the wrong sign. This is due to the fact that 
(13) has been derived from an antiferromagnetic Heisenberg hamiltonian. Although this 
wrong sign apparently produces an instability in the ferromagnetic spin wave spectrum, this 
instability is not significant. Recall that the ferromagnetic phases in doped manganites are 
due to the interaction with the charge carriers. Hence any reliable estimate of the spin 
wave dispersion relation parameters in the ferromagnetic phase must also take into account 
the interaction with the charge carriers. We shall do so later on. Such kind of (fictitious) 
instabilities also occur in the canted phases although they are not so immediately spotted 
from the Lagrangian ( [4.6| ). 



4.1.2 Antiferromagnetic configuration 

In the antiferromagnetic configuration S ^ and VL —>■ 2M, and after performing the 
corresponding rotation, 1 — > 3 — > 2, as explained in section 2, we obtain the following 
effective lagrangian: 



GJAFCi^ 



a, a,' 



(4.^ 
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In this case, since the ground state configuration is supported by the Heisenberg hamil- 
tonian, the spin waves obtained from it are stable. 



4.2 Spin waves dispersion relations: contributions from charge 
carriers 

We have now at our disposal suitable low energy effective lagrangians which describe spin 
waves in doped manganites in the canted, ferromagnetic and antiferromagnetic phases. They 
are given by the pure spin wave terms above together with the terms of interaction with the 
charge carriers (|3.16| ). More precisely, the effective lagrangian for the canted phases can be 
obtained from ( |3.16| ), ( p.l7[ ) and ([4.6|), for the ferromagnetic phase from ( |3.16| ), ( p.2(]| ) and 
( [4.7| ), and for the antiferromagnetic phase from ( |3.16| ), ( p.24| ) and 



In order to obtain a reliable evaluation of the parameters in the spin waves dispersion 
relations we have to take into account the interaction with the charge carriers in the spin 
waves two point Green's functions. This can be easily achieved from a further (this time 
non-local) effective lagrangian which is obtained by integrating out the charge carriers and 
keeping only the contributions up to two spin wave fields. 

By integrating out the fermionic fields in the equation ( p.K ) we obtain the following 
contributions to the effective lagrangian: 

5g = -iTrlogiL^+Wn) = -iTr log L^-iTr{L:[^Wu) + '-Tr{L^'WuL^'Wu)+- ■ -, (4.9) 

where Tr stands for the trace over the space-time indices as well as the matrix indices. We 
have expanded the logarithm up to second order. 

The first term in ( [4.9| ) gives rise to an effective potential for the ground state configuration 
which, together with the static antiferromagnetic Heisenberg term, produces the rich phase 
diagram for doped manganites presented in [^. The following two terms in the expansion 



are responsible for the appearance of terms with at least two spin waves in the effective 



lagrangian, as can be see from ( 3.1 



Even though in order to obtain the relevant contributions to the effective lagrangian it 
is enough to consider the interaction up to two spin waves in ( |3.18|) , interesting general 



characteristics will arise if, instead, we use the SU{2) invariant expression (|3.17|) for the 



interaction. In this way we are going to obtain not only an explicitly invariant effective 
lagrangian under SU{2), but also the non-local structure which arises from the absence of 
gap in the fermionic spectrum of excitation. 

We begin with the calculation of the second term in (^4.9| ), i.e. = —iTr{L^^Wii). 

In this calculation a closed loop integral, representing the density of carriers, appears, 

X , r°° duj dk 
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X is the doping and the volume of the unit cell, tr represents the matrix trace, the 
space-time trace has already been taken into account in the integration over the momentum 
q = (a;, k). Lj~^(g) is the Fourier transform of the inverse of the operator Li given in ( p.l3| a). 
The convergence factor e*'^'' [r] — > 0"^) is introduced to pick up the correct order of the fields 



in the calculation of closed loops of one point Green's functions [^] . 
The contribution of these terms to the effective lagrangian is given by 

1, 



(4.11) 



2a3 2 ° 2ma3 

where we have dropped terms which contribute with a total derivative. 

Whereas all the contributions in ( ^.11| ) are local, because the loop integral is closed, the 
contribution from ^^J^^ = ^Tr{L^^WnL^^Wn) is going to contain non-local terms due to 
the presence of the so-called vacuum polarization tensor 

nS^'^(p) = -'Jj^Ap + lyLiaiP + Q) <l'LY,\q), (4.12) 

where a, 6 = +,— represent the diagonal components of the operator Li given in ( p.l3| a). 
i,j = 1,2,3 represent the spatial components of the momentum, while i,j = means 
the absence of the corresponding momentum component. The properties of this tensor are 
displayed in the Appendix. 

Taking into account the symmetry properties of the vacuum polarization tensor the con- 
tribution to the effective lagrangian reads 



0(2,2) 
^eff 



dudw 



dp 



-ip{u—w) 



(4.13) 



■ cos 



1 



m 



b-(u)b+(w) + bt(u)bj(w) + bt(u)b+(w) + br(u)bj(w) 



where summation convention over repeated indices has been used, and as in the previous case 
terms contributing with a total derivative to the effective lagrangian have been dropped. It is 
easy to see that this part contributes with non-local terms as long as the vacuum polarization 
tensor has a non-constant behavior in the energy-momentum vector = {u, p). One of the 
most interesting terms with these characteristics is h^{u)h\{w), which mixes time and spatial 
derivatives. 

The leading contribution to the effective lagrangian is given by keeping in ([4.1 3|) second 
order terms in derivatives (or momentum). This corresponds to consider the zero energy 
and momentum limit of the vacuum polarization tensor, i.e. n^*^''^(0). It is also convenient 
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to choose the basis ([A.4|) , which has its third component parallel to the momentum, p, in 
order to simplify the calculation. In this basis, and using the relations given in (|A.8|) and 
( |A.9| ), the action ( [4 .131 ) reads 



(2,2) 



J dudw J 



e-'P("-"')|-cos2 - n(°'°)fO) 



e X 



- cos - 



2 ma;^ 



(e;3)^'(«) 



+ 77 COS 



e 



(4.14) 



+ {e\^)ht{u)) {el)b^{w)) + (ej„)6r(«)) {e\^^b.{w 



Hence, with the aid of ( p. 12 ) and ( p. 11 ) we can expand this expression up to two spin 

), i.e. 



waves. At this order only one spin wave must be considered in the expansion of (|2.11 



6^ ~ <9^7r", which means that they are proportional to the energy- momentum, and since the 
vectors e(i) and e(2) are perpendicular to the momentum, e^^-jfo" ~ ^\a)'^iT^'^ ~ ^e^^-^pV* = 0, 
(a = 1, 2), they do not contribute at this order. 

In addition to this, it is very interesting to notice how the terms in the second line, which 
would contribute with time derivatives for vr^f-u), cancel at this order, 



IpI P 



(4.15) 



The cancellation of these terms is very important, since as it can be seen from (|A.8| ) the 
tensor n^°'°)(0) contains an imaginary part, which would produce the spontaneous decay of 
the spin wave 7r^('u) into fermionic excitations. 

The final result for the effective lagrangian up to two spin waves, after using a similar 
procedure to (|4.15|) for spatial derivatives, turns out to be 

JafM"^ 



£(2,2) ( 



U 



2mt/2 



n-i__ sin^ - 



+ 



X o 



cos - diTf diTl 



(4.16) 



where n+_ is given in ( |A. 12| ) and ( |A.14| ) for one band canted {CCl) and two band canted 
{CC2) phases respectively. 
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4.2.1 Ferromagnetic configuration 



The ferromagnetic limit is again very easily taken from the canted results when 6^0 
(equivalently we could use ( p.20D ). In addition to that, we notice from ( p.22| ) that the 
interaction already contains at least two spin waves, thus in order to calculate the dispersion 
relation we only have to consider the second term in the expansion of the logarithm 
The effective lagrangian at this order is 



4.2.2 Antiferromagnetic configuration 



X 



2ma^ 



a, a; 



(4.17) 



In the extreme low energy and momentum limit we are interested in, there will be a 
contribution for the antiferromagnetic state in the conducting phase only. As in the previous 
cases the antiferromagnetic limit of the canted expression must be taken carefully. In this 
case, as in the ferromagnetic one, the interaction already contains at least two spin waves 
and it is enough to consider the second term in (|4.9| ). Since the antiferromagnetic state 
corresponds to a 2-fold band, Li is degenerated, the term proportional to in ( |3.24| ) will 
cancel, which prevents a term with a single time derivative from appearing in the effective 
lagrangian as it should be. The final result reads 



7^ a- 



(4.18) 



4.3 Spin waves dispersion relations: final results 

Finally by summing all the contributions C = C^^^ + C^'^'^^ + C^'^''^\ given in (|4.6| ), ( [4.11|) 
and ( [4.1(j| ) respectively, we are in the position to write the effective lagrangian for the spin 
waves up to second order in derivatives, and up to two spin wave fields. After expanding the 
first two contributions in spin wave fields and taking into account the expression (|2.17|) we 
obtain for the parameters of the spin waves the following results: 
2M + x 



B 



C 



fl 



,3/2 



hA ( Qtt'^x 



2m ISvrM (2M + x)y 

1 z^/^ 1 
2m 157T^A (2M + x)y 



[l + n+_](l-y2 



.3/2 



^3/2 



2mz ISttM M2 



2 + 



5 A I Gtt^x' 



x3/2 



(4.19) 
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where y = cos^/2 is a measure of the canting angle 6 and x is the doping. n+_ also depends 
on X. We have used the expression of the parameter A in ( |A.llj ) and 2m ~ z/a^t. In the 
case of the one band canted phase (CCl) n+_ is given by ( [A.12| ), and all the expressions can 
be written explicitly in terms of the doping using ( A. 11 ). In the case of the two band canted 
phases (CC2) n+_ is given by (|A.14|) , but it is impossible to write all the above expressions 
explicitly in terms of the doping only; we need also the canting angle, y, which depends 
implicitly on the doping. For a given value of the doping we can obtain the corresponding 
value of y by solving the equations (|A.13|) . 

As it was described in section 2 a Bogolyubov transformation must be carried out in 
order to diagonalize the lagrangian and obtain the physical fields, which have a mass given 
by l/2m' = - in (|A|), 



1 



,3/2 



2m' 2m I^tt^A (2M + x)y \ 



5A ( Gtt^x 



y3/2 



/ Gtt^x' 



r3/2 



-22/2_2n+_(l-y2). 

(4.20) 

The expressions (|4.19|) and ( [4.20|) for the velocity and the mass of the spin waves are 
used in the plots of Fig. 1. 



4.3.1 Ferromagnetic spin waves 



The parameters for the ferromagnetic spin waves must fit (p.24|) , and are obtained from 
the sum of (g^) and ( ^T^ ). 

2M + X 



f 

1 



a-" 

1 ^3/2 



2m' 



-2 + 



5 A ( evr^x' 



,3/2 



(4.21) 



2m ISvrM (2M + x) 

which corresponds to the limit of the canted parameters ( [4.19|) taking into account that in 
the ferromagnetic limit l/2m' = B. This is plotted in Fig. 1 in the case of a ferromagnetic 
phase. 

In this limit configuration it is particularly easy to see the effect of the charge carriers in 
the behavior of the spin waves. Since the interaction between the core spins is antiferromag- 
netic the mass derived from this interaction is negative (producing an unstable spin wave). 
However, the contribution from the fermionic sector compensates this sign and stabilizes 
the spin wave. It is also worth noticing that in this limit from the canted phase one spin 
wave disappears, namely, t^^{u). Due to the existence of the remaining f/(l) symmetry in 
the ferromagnetic ground state this fluctuation does not modify the ground state any more. 
We can see from ( |4.19| ) that the velocity of the linear branch goes to zero smoothly in the 
transition to the ferromagnetic phase. This answers the question asked at the end of the 
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first paper in |3^: "How do the spin dynamics evolve from tlie double-brancli state reported 
liere into a state witli only one branch characteristic of a ferromagnetic metallic phase?" . 



4.3.2 Antiferromagnetic spin waves 

In this case the parameters for the effective lagrangian are obtained by summing the 
contributions of 



f 



and ([4.18|) , and are given by 

z 



^3/2 



2mz Ibn^A 




(4.22) 



which coincides with the limit of the parameters for the antiferromagnetic spin wave in the 
canted phases, 7t^{u). In the insulating phase we should take x = 0. This velocity is plotted 
in Fig. 1 for the antiferromagnetic phase. 

As in the previous case we should notice the lost of one spin wave field when we carry 
out the limit towards the antiferromagnetic ground state, 7r^(u) (after the rotation of the 
reference system, 1 — 3 — > 2, it becomes n^{u)). The combination of the canted phase 
parameters f^{B + C) with y ^ 0. It is interesting to notice that the other remaining 
field, which used to be part of a ferromagnetic mode, with a quadratic dispersion relation, in 
the canted phase, becomes part of the antiferromagnetic spin wave with a linear dispersion 
relation. The remaining combination of the canted phase f^{B — C) — » when y ^ 0. 
This answers the question at the end of the previous subsection in the antiferromagnetic 
limit. 



5 Disentangling Canted Phases from Phase Separation 
Regions 

Recently controversial results have appeared in the literature regarding the existence of 
canted phases in doped manganites, in particular concerning their stability against phase 
separation ^ ^ We showed in that canted phases not only exist but they are 
also thermodynamically stable. We presented there a phase diagram where, in addition to 
stable canted phases, phase separation regions appear. 

The phase diagram presents the following phases: antiferromagnetic insulating (AFI), 
antiferromagnetic conducting with two bands {AFC2), canted conducting with two bands 
(CC2), canted conducting with one band (CCl), ferromagnetic conducting with one band 
(FCl) and four phase separation regions between the FCl phase and the remaining, i.e. 
PSl {AFI - FCl), PS2 {AFC2 - FCl), PS3 [CC2 - FCl) and P^4 [CCl - FCl). 
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The question of which regions the system passes through when going from the antiferro- 
magnetic insulating phase to the ferromagnetic conducting one upon increasing the doping 
in the actual materials could not be satisfactorily solved there, since the answer depends 
critically on the values of the parameters of the model. For reasonable values of these pa- 
rameters, various possibilities are allowed. We have chosen five values of the parameter A, 
in (|A.11| ), for which the sequence of phases is the following: 



A = 


2.20 


AFI - 


PSl - 


FCl 






A = 


1.75 


AFI- 


AFC2 


-PS2 


-FCl 




A = 


1.40 


AFI- 


AFC2 


-CC2 


-PS3 


-FCl 


A = 


1.00 


AFI - 


AFC2 


-CC2 


-CCl 


- PSA 


A = 


0.80 


AFI - 


AFC2 


-CC2 


-CCl 


-FCl 



In order to establish differences between the canted phases and the phase separation 
region we must make a guess on how these phase separation regions look like, since our 
model does not describe these non-homogeneous regions of the phase diagram. Even though 
these may be very rich regions, with many different structures in them, as charge ordering, 
stripes, orbital ordering or polaronic excitations pR O, H^, 03, EO, pOl E^, we shall assume 



that the main structure is the coexistence of two macroscopic domains corresponding to the 
phases at the border of the phase separation region, and that the interphase will not disturb 
qualitatively the properties of each of them. With this assumption in mind we have plotted 
in Fig. 1 the dependence on the doping for the velocity and the mass of the spin waves for 
each of the values of the parameter A given above. 

The dotted lines corresponds to the values of the velocity and the mass in the phase 
separation region, for the first four values of the parameter A. They are constant values, 
because they are given by the value of the corresponding phase in the border of the phase 
separation region. Consider, for example, the first case, where an AFI and a FCl domain 
coexist in the phase separation region. The doping is an extensive magnitude, and even 
though it reduces globally over the system, the density of carriers remains constant in the 
FCl domain, since this one reduces as the doping (considered globally) decreases. 

In the third and fourth cases, where canted domains coexist with the ferromagnetic 
domain, we can observe two different values for the masses of the spin waves in the phase 
separation region, as well as one velocity. 

Let us concentrate in the first versus the last case, i.e., the phase separation region with 
antiferromagnetic insulating and ferromagnetic conducting {AFI — FCl) domains versus 
the canted phases [^, ^ ^]. Since the differentiation between these two structures 

seems to be an experimental challenge, we describe below a few distinct properties of the 
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V2 



m/m' 




Figure 1. The dependence of the velocities and the masses with the doping for five differ- 
ent values of the parameter A (~ t/JAp)- = {Ihir'^ A)2mzv^ /Qz^/'^{JAFa^M'^) and in/m! = 
{Ihir'^ A)m/ z^^'^m' . The horizontal dotted lines correspond to the phase separation regions, and the 
vertical dotted lines correspond to the phase transitions. 
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spin waves which may help to differentiate between a canted phase and a phase separation 
region consisting of ferromagnetic and antiferromagnetic domains: 

(i) First of all, in the canted phases dispersion relations we observe one ferromagnetic 
branch and one antiferromagnetic branch, whereas in the phase separation case 
{F — AF) we should observe one ferromagnetic, but two antiferromagnetic branches. 

(ii) Second, the antiferromagnetic branches present a further dramatic property: its be- 
havior in the presence of a magnetic field along the staggered magnetization. Whereas 
in the antiferromagnetic case the two branches will be splitted, in the canted phase the 
single linear branch will not be even shifted by the presence of such a magnetic field. 

(iii) Finally, we have presented in the previous section and in Fig. 1 the different behavior of 
the dispersion relation parameters with the doping, A = 2.20 for the phase separation 
region and A = 0.80 for the canted phases. 

These three characteristics, in particular, the first and the second one which are of rather 
general nature (model independent), should allow to experimentally differentiate the regions 
of the phase diagram where ferromagnetic and antiferromagnetic phases coexist from those 
where real canted phases exist. 



6 Conclusions 

We have presented a complete study of the spin waves in canted phases. We have ex- 
ploited the spontaneous symmetry breaking pattern SU{2) ^ 1 to construct an effective 
lagrangian for low energy and momentum spin waves in canted phases at next to leading 
order. For simplicity, we have chosen a cubic lattice, but any other lattice can be treated 
within the same formalism. The lagrangian at leading order depends on five parameters and 
to next to leading on nine. The leading lagrangian yields two spin wave modes, one with 
a quadratic and one with a linear dispersion relation. The leading effective lagrangians for 
ferromagnetic and antiferromagnetic ground states were also considered as limit cases of the 
canted configuration. These depend on two parameters only. 

Since the canted phases appear in doped manganites, and are associated with conducting 
properties of these materials we have also presented interaction lagrangians of spin waves 
with charge carriers. Whereas the lagrangian for spin waves alone is of general nature 
(model independent), the interaction with charge carriers depends on microscopic features 
of the material, in particular on the number of conducting bands available. We have chosen a 
simple case (two band) which is inspired in a realistic model introduced in for the study 



of the phase diagram of doped manganites. Again we have derived the interaction lagrangian 
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for ferromagnetic and antiferromagnetic phases as the extreme cases of the interaction with 
canted spin waves. 

We have apphed our results to the study of the spin waves in doped manganites using the 
continuum double exchange model |^ . We obtained the explicit dependence on the doping 
and the canting angle for /^r, l/2m', and v, which determine the dispersion relations for 
the spin waves, and we have plotted them in Fig. 1 for several parameters. 

Finally, we have proposed three ways to tell experimentally apart canted phases from 
phase separation regions (coexistence of ferromagnetic and antiferromagnetic phases) by 
looking at suitable properties of spin waves. The results above may also be useful for a more 
refined study of the phase diagram of doped manganites. In particular the leading quantum 
corrections to the classical spin dynamics in the low energy region are due to spin waves. It 
would be very interesting to elucidate the effect of these corrections in the phase diagram. 
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Appendix A: Vacuum Polarization Tensor 

The vacuum polarization tensor, as was defined in subsection 4.2 reads 



where 



L^a\'l) = ^ , fi± = -^^^l + y±27cos^-/i. (A.2) 

a; - Qa + if]^ 

2m 

The symmetry properties of this tensor under the change of sign of the energy and momentum 
are given by 

nlf(^,-p) = (-!)'-+''« nl^^V,p) (A.3) 
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In order to simplify the calculation of the integrals we chose a reference system with its third 



component parallel to the external momentum in Il^^^''\u,p), namely, 



6(1) 



(pk)p - k 
|k X p| 



3(2) 



k X p 

|k X p| 



6(3) = P, 



(A.4) 



where k is a unit vector in the third crystallographic direction. 

We will denote the components of the vectors and tensors in the new basis with Greek 



indices, instead of Latin ones, such that they verify 



Wi e*^^-), and consequently. 



(A.5) 



Therefore the symmetry properties under the change of sign of the energy and momentum 
in the new basis become 



n 
n 



ab ' 
ab ' 



n 



(",/3) , 
ab I 



i^, -P) 
-P) 



(_i)5o,+5o/3 n[f'"V>p) 



(A.6) 



In this basis, after some straightforward algebra, following [Q, it is easy to see that all the 
dependence on the components and of the integrand in ( |A. 1| ) comes from the explicit 
dependence [p + q)" and q^, and therefore it is verified that 



(A.7) 



where the repeated indices are not summed. 



Once the symmetry properties of the vacuum polarization tensor have been considered 
we will address its calculation. As it was stated in the subsection 4.2 it is enough to calculate 
the limit — > of the tensor, i.e. n^*^"'''(0), to obtain the leading contribution to the effective 
lagrangian. Since the external energy, u, and momentum, p, injected in the fermionic loop, 
are related, this limit must be taken carefully. If they correspond to a 7r'^(a;) spin wave their 



relation is linear, u 
components are given by 



p|, while the relation for 7r^(a;) is quadratic, z/ ~ p^. The tensor 



nr(o) 

nif)(o) 

niy)(o) 



47]-2 



-2 + Xa log 











-in\xa\9il - \Xa\) 



Xn 



ka\p\ 



^,o{-n^) + hl{i-xl) nir(o) 



(A.8) 
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nf_:")(o) 



+ nr(0) 

m kle{-Q+) - k^_0{-Q_) 



157r2 



ttiJafM'^ 
t 



(a = 1,2,3), 



where ka = ^ —2m£la represents the Fermi momentum in each band a = +, — , and Xa 
is the relation between the spin wave velocity and the Fermi velocity in each band. The 
step functions 6{—Qa) ensures that only the bands which are bellow the chemical potential 
contribute to the result. In the expressions above the summation convention was not used. 
An interesting result arises from considering the summation over n^^'^-'(O), namely. 



—m- 



{2m 



1 3/2 



E(-^J^/'^(-^.) + ff^)'Enir(o) 

a \ IPI / a 



mx 



(A.9) 



as can be verified from an explicit calculation of the integral which gives the doping x in 

(EE)- 

The results obtained until now are exact (in the limit of 0). In order to obtain 

further analytic results for II^l"'' (0) in ( [A.8|) we will use the value of Q± to leading order in 
7, already obtained in p2| , 

fi± = -%o±Z/), (A. 10) 

where i/q is a measure of the chemical potential. 

In the CCl phase we obtained in p2[ the following values for the canting angle, y, and 
the doping, x: 



A 



.3/2 



t 



157r2 JApa^M^ 



(A.ll) 



and therefore 



CCl : 



a^A kl 



.3/2 kl-kl 



4 / Gir'^x 



2/3 



5 V ^3/2 



(A.12) 



The calculation is a little bit more complicated in the two band case, CC2, where we 
obtained in P2[ for y and x the following expressions: 



y 



-A 



{yo + y?^' - {yo - yf' 



\A{y' + ^yl) = iyo + yf' + iyo-yf' 



(A.13) 



29 



which yields after ehminating the chemical potential measure, yo, 



CC2 : n , 



j3/2 kl - 



1 4 



4 f67c^x\ ^ A f 67r^x\ y-.o 5A f 6'K^x\ , 



v/35\ 5A \ ^3/2 J ^ ' 2 



where the limit y ^ can be taken smoothly, and the previous result is obtained by taking 

1 /2 

into account that Ay^' =8/15 when y = 0. 
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